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Abstract 



(N 

■ We study the symmetry behavior of the Gross-Neveu model in three and 

o: 

two dimensions with random chemical potential. This is equivalent to a four- 
fermion model with charge conjugation symmetry as well as Z2 chiral sym- 
metry. At high temperature the Z2 chiral symmetry is always restored. In 
three dimensions the initially broken charge conjugation symmetry is not re- 
stored at high temperature, irrespective of the value of the disorder strength. 
In two dimensions and at zero temperature the charge conjugation symme- 
try undergoes a quantum phase transition from a symmetric state (for weak 
disorder) to a broken state (for strong disorder) as the disorder strength is 
varied. For any given value of disorder strength, the high-temperature behav- 
ior of the charge conjugation symmetry is the same as its zero-temperature 
behavior. Therefore, in two dimensions and for strong disorder strength the 
charge conjugation symmetry is not restored at high temperature. 
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I. INTRODUCTION 



Intuitively, when heated, a system with initially broken symmetry will recover its symme- 
try because thermal fluctuations are able to overcome potential barriers. But a counterex- 
ample was noticed by Weinberg [[[]]: For the four-dimensional 0(N)xO(N) scalar 4 model, 
he showed that the system can remain in the broken phase even at sufficiently high temper- 
ature. This phenomenon is called inverse symmetry breaking or symmetry nonrestoration 
(SNR), depending on whether the system was in a symmetric or a broken phase at zero 
temperature. 

Since Weinberg's observation, SNR has been a subject of academic curiosity or a can- 
didate way out of cosmological problems caused by topological defects like monopoles and 
domain walls (see Ref. [@] for review). According to Bajc's classification ||, there are three 
classes of SNR mechanisms in field theory: (i) a prototype case like the two-scalar model 
IHH, (ii) flat directions in supersymmetric theories |4j], and (iii) large charge density (or 
chemical potential) |5]||. Here we restrict ourselves to class (iii). 

If a large enough charge can not be stored in thermally excited modes at high tempera- 
ture, it must reside in the vacuum, and this is a sign of SNR. In field theory, a scalar field 
(order-parameter field) gets a positive mass term by thermal effects, but a negative one by 
the effects due to the chemical potential. For a fixed charge (i.e., in the canonical formalism), 
the chemical potential is temperature dependent. In this case, if the effect of the chemical 
potential on the mass exceeds the thermal effects at sufficiently high temperature, the scalar 
field acquires a nonzero vacuum expectation value (i.e., SNR) ||. 

However, in an open system of the model which does not belong to the class (i) or (ii), the 
symmetry may always be restored at high temperature. In the grand canonical formalism, 
the chemical potential and the temperature are independent parameters and so the thermal 
effect on the mass always surpasses the effect of the chemical potential at sufficiently high 
temperature, for a fixed chemical potential. For example, consider the Gross-Neveu (GN) 
model ]7^|§ with chiral symmetry. At finite chemical potential the initially broken chiral 
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symmetry is always restored at high temperature |10 . 

In order to find a new kind of SNR in four-fermion models, we will extend the GN model 
at finite chemical potential to a disordered model with random chemical potential. Recently, 
disordered nonrelativistic Dirac fermions in two spatial dimensions have been studied in 
relation to the integer quantum Hall transition [[□]]. Pure fermions exhibit such a transition 
as the value of the mass is tuned through zero, but its universality class is different from the 
one observed in actual experiments. Usually three types of (static) disorder are considered 
for a more realistic model: random gauge potential, random chemical potential and random 
mass. 

Motivated by the SNR mechanism (iii), we introduce the (relativistic) GN model with 
random chemical potential in Sec. II. If the chemical potential has a Gaussian distribution 
at each site, our model is equivalent to the four-fermion model with two kinds of four- 
fermion interaction, (\I/\]/) 2 and (\I/7o\I/) 2 (see Eq.(3)), and has charge conjugation symmetry 
in addition to the Z 2 chiral symmetry. In Sec. Ill we examine the behavior of these sym- 
metries as the temperature or disorder strength is varied using the 1/N expansion in three 
and two dimensions. While Z 2 chiral symmetry is always restored at high temperature, the 
charge conjugation symmetry exhibits SNR. In addition, we check the validity of the mean 
field approximation (the leading approximation in the 1/N expansion) in two dimensions. 
In Sec. IV the fundamental origin of SNR for charge conjugation symmetry is discussed 
conceptually. Our conclusions are presented in Sec. V. 

II. GROSS-NEVEU MODEL WITH RANDOM CHEMICAL POTENTIAL 

The Euclidean Lagrangian of the GN model at finite chemical potential \i is given by 

£ = 1 |r( /9 + yU7o )M/_^-(M) 2 , (1) 

where g 2 {> 0) is the coupling constant of the four-fermion interaction (\i/\I/) 2 and N is 
the number of flavors of the Dirac fermion The 7 matrices are 4x4 and hermitian. 
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Let us consider the system under the influence of a random chemical potential p{x) with 
the Gaussian distribution exp(— / d d x^2p 2 ) at each site where R 2 (> 0) is the strength of 
disorder and d the dimension of the Euclidean space. The Gaussian noise is characterized 
by correlation functions 

(p(x)) = 0, (p(x)p(x')) = ^5 d (x-x f ). (2) 

After integrating out the random chemical potential, our model is equivalent to the 
four-fermion model 

C = ¥ /M - ^ [^(M) 2 + i? 2 (% ^) 2 ] , (3) 

with the Z 2 chiral symmetry — > 75^, — > —^75} and the charge conjugation symmetry 
_^ cW, * -> -^ r Ct}. Here the matrix C satisfies C7+C = 1, C^C = - 7 J. Under 

charge conjugation, and ^70^ transform to and — ^70^ respectively. Hence, the 

Lagrangian Eq.(l) with definite chemical potential p does not possess the charge conjugation 

symmetry (i.e., fermion-antifermion symmetry). Note that in Eq.(3) the chemical potential 

term does not appear explicitly. 

We will study this model by the leading approximation of the 1/N expansion in three 

and two dimensions. To easily incorporate the 1/N expansion, let us rewrite the Lagrangian 

Eq.(3) by introducing scalar auxiliary fields a(x) and p(x): 

£ = ^(f) + a + p 7 o)* + ^- 2 a 2 + ^p 2 . (4) 

The random chemical potential p(x) plays the role of a scalar auxiliary field. The Z2 chiral 
symmetry and the charge conjugation symmetry are now expressed as — > 75^, ^ — > 
-*7 5 , a -> -a} and -> C^ T , * -> -^C 1 ", p -> -p}, respectively. 

III. THE BEHAVIOR OF Z 2 CHIRAL SYMMETRY AND CHARGE 
CONJUGATION SYMMETRY AT ZERO AND HIGH TEMPERATURE 

For finite-temperature field theory we adopt the imaginary-time formalism. At inverse 
temperature (3{= T~ l ), the fermion fields are antiperiodic on R d ~ l x [0,(3], while the scalar 
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auxiliary fields are periodic. Let us introduce the notation: f^ T ' = T Y£?=-oo I (2 7 r) d - i ' 



r(°) = 

Jo — 



/ jljSx- Integrating out the fermion fields in the partition function for Eq.(4) we obtain the 
effective action for the auxiliary fields a and p. In order to investigate the vacuum structure 
we need to find the finite-temperature effective potential Vt((t, p) by taking a and p as 
constant fields: To leading order in the 1/N expansion, 

,(T) 



Vt((t, p) = 

N 2g 2 2R 2 



/v ) 

2 / ln[(po - ipf + p 2 + a 2 

Jv 



(5) 



where p$ = (2n + 1)^/(3 = uo n (n = integer) at nonzero temperature. Note that the effect of 
the chemical potential p is to shift the energy by —ip. 

To evaluate the integration in Eq.(5), we need some mathematical formulae. By the 
standard method of contour integration fi~3|| , 
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uj n - ip 
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(u n - ip) 2 + cr 2 



ipT E 



l + e /3(M+H) 1 + e 0(M-H) 

< ~ (* 2 ~ P 2 ) 
(u 2 + a 2 - p 2 ) 2 + (2pu n ) 2 

sinh(/3p) 



(6) 



(?) 



cosh(/5a) + cosh(/3p) 

When the GN model is studied in the canonical formalism (i.e., with a fixed charge 
similar calculations appear with imaginary chemical potential. In this case a regulating 
factor of the form e lu)nT is needed in evaluating the summation in Eq.(7) and ensures a finite 
result in the limit r — > after the Matsubara sum has been performed. By using Eqs.(6) 
and (7), we obtain 



T E MK - W? + ° 2 }=T [In 2 + ln(cosh(/5a) + cosh(/3p))] , 



(8) 



where the (^-function regularization was used to determine the field-independent constant. 
At zero temperature, these formulae reduce to 
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2ir (p — ip) 2 + a 2 
dpo Po - ip 



fl(H-H) 
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2n (po - ip) 2 + a 2 = 2 Sgn(/?)W ~ ^ 



[ ln[(po - ip) 2 + a 2 ] = max(|(r|, \p\). 
J 2ir 



(9) 
(10) 
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Using Eqs.(6), (7) and E a = VP 2 + a 2 , we have 



^- ^ = ^ - 2a 



da \ N J g 2 J (2ir) d ~ 1 E a 



d fV T \ a f d^p 1 



1 - 
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1 + e P{E a +p) I _|_ e P{E a -p) 



(12) 



dp\N) ~ R 2 smh ^P> J ( 27r )d-i cosh(/3K) + cosh(/3p) " 1 dJ 

To renormalize the effective potential Vr(a, p) let us consider the GN model at zero 
temperature and in the absence of the random chemical potential because the effects of 
temperature and chemical potential do not change the ultraviolet behavior. Define 1/G 2 = 
1/g 2 - 1/g 2 with 1/g 2 = 4 / In 2 < d < 4, the GN model is in the broken phase of 

the Z 2 chiral symmetry for negative G 2 , corresponding to strong coupling (g 2 > g 2 ) in the 
cutoff regularization, while it is in the symmetric phase for G 2 > 0, corresponding to weak 
coupling (0 < g 2 < g 2 ). In particular, in two dimensions, the Z 2 chiral symmetry of the GN 
model must be broken no matter how we choose the coupling g 2 ||. In the broken phase, 

1 = 4 1 (14) 

g 2 J (27r)V + M 2 ' 1 J 

where M = |(er)|(> 0) is the dynamically generated fermion mass at zero temperature. 

From now on, we will adopt dimensional regularization, where G 2 is equal to the regu- 
larized g 2 . 

A. Three dimensions 

In this case, renormalization is not needed to the leading order of the 1/N expansion (in 
dimensional regularization). By making use of Eq.(ll), we can find the zero-temperature 
effective potential Vo(a, p) directly: 

V (a,p) a 2 p 2 1 r 3 2 , , .."I 

— N~ = 2G 2 2R 2 ~ 6^ ^ maX ^ ~ max (l a l' \P\)\ ' ( 15 ) 

where 1/G 2 = —M/tt for broken Z 2 chiral symmetry. The gap equations have four kinds 



of solution (\a\,\p\): (i) (0,0), (ii) (M,0), (hi) (0, 2n/R 2 ) and (iv) (y M(M — 2n/R 2 ), M). 
The solution (iv) exists only for M > 2tt/R 2 and corresponds to saddle points. Fig. 1 shows 
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the zero-temperature effective potential as a function of a/M and p/M, for broken Z 2 chiral 
symmetry, (p) = is metastable, irrespective of the values of G 2 and R 2 . For \p\ > \a\, 
however, Vo(cr, p) is unbounded from below due to the — |p| 3 /(67r) term, which indicates 
breaking of the charge conjugation symmetry. This result stems from the fact that the term 
— |p| 3 /(67r) arising from quantum effects surpasses the effect of the probability distribution 
(p 2 /(2R 2 )) for large \p\. 

At finite temperature, using Eqs.(6), (7) and dimensional regularization, we obtain 



(T) 1 1 

u n - ip) 2 + El = ~4^f3 
( T ) u n - ip _ i sgn(p) 
(u n - ip) 2 + El ~ 47T/5 2 



C j^wtei = ~h ^ 2e " w coshm + e "" w )l ' (16) 

« ln H + U 2 (-e^) - Li 2 (-e^H„))| . (17) 



Here the polylogarithm Li u (z) is defined (for v > 0) as Li u (z) = J2t=i TP ( see ^ e ^- E3 ^ or 
useful properties). From these formulae, the finite-temperature effective potential Vr(cr, p) 
is given by 

Yl^lA = ^ + ^_^! + A_ [Li,(-e^)) + LU-e^-ri) 
N 2G 2 + 2R 2 3tt + tt/33 L M } + 31 } 

-Pa {Li 2 (-e /3 ^) + Li 2 (-e^))}] , (18) 

up to a field- independent constant. At sufficiently high temperature, 

V T (a, p) /In 2 



N \ 7T 



7V-P 2 )- (19) 



While the initially broken Z2 chiral symmetry is restored at high temperature, charge con- 
jugation symmetry is not. Hence our model exhibits nonrestoration of charge conjugation 
symmetry irrespective of the values of G 2 and R 2 . We may interpret this phenomenon as an 
inverse symmetry breaking because (p) = is metastable at zero temperature. Intuitively, 
SNR is related to the tachyon-like behavior of the random chemical potential (see Eqs.(19) 
and (24)). In the quantum correction term of Eq.(5) the chemical potential acts as a nega- 
tive mass term (— p 2 ) contrary to the usual positive mass term (cr 2 ). From a different point 
of view, we will discuss the origin of SNR conceptually in Sec. IV. 
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B. Two dimensions 



For dimensional regularization, we work in 2 + e dimensions. In terms of the fermion 
mass M, the zero-temperature effective potential is given by 

max(|a|, \p\) + y / max 2 (|(r|, |p|) - a 2 



V (a,p) _ a 2 



N 



2tt 



-1 + 2 In 



M 



2R 2 



— V max 2 (|ct|, |p|) -a 2 , 

71 v 



(20) 



where Eqs.(ll) and (14) were used. Unlike in three dimensions, for large \p\ the effect 
of the probability distribution (p 2 /(2R 2 )) is comparable to the last term —p 2 jix) in 
Eq.(20) arising from quantum effects. The charge conjugation symmetry can be controlled 
by the strength of disorder R 2 . The system is in the symmetric state for < R 2 < tt/2, 
while in the broken state for R 2 > tt/2. Fermions and antifermions are equally proba- 
ble in the symmetric state ((p) = 0), but only fermions (or antifermions) are allowed in 
the broken state ((p) = ±oo). Our system suffers from a quantum phase transition at 
R 2 = 7r/2(= R 2 C ). The gap equations have solutions (|cr|, \p\): (i) (0,0), (ii) (M, 0) and (hi) 



{^{2R 2 - n)/(2R 2 + tt)M, 2R 2 M/(2R 2 + tt)) for R 2 ^ R 2 C , and (i) (0,V|p|) and (ii) (M, 0) 
for R 2 = R 2 . The solution (iii) exists only for R 2 > R 2 and corresponds to saddle points. 
Fig. 2 shows the zero-temperature effective potential as a function of a/M and p/M in (a) 
the symmetric and (b) the broken phase for the charge conjugation symmetry. 

To examine the high-temperature (f3 — > 0) behavior, let us introduce dimensionless 
quantities: Vr = P 2 Vr, o = (3a, p = ftp, M = f3M. We want to expand the finite- 
temperature effective potential in a and p at high temperature (i.e., for small a and p). In 
terms of the dimensionless quantities, Eqs.(12) and (13) are reduced to 

1/1 1 



d_ Vr 
da [ iV 



d (Vt 



2a 

7T 

2a 

7T 



In 



\a\ 



M 



dx- 



+ 



In 



^Jx 2 + a 2 V 1 + e Vx ' 2+a2 +P 1 + e^ x2+a ' 2 -P 

(X)- 7 + o(^)} + {^) + o(. 2) ) p2l - 



0(P 4 



(21) 



1 2sinh(p) f°° 

'r 2 ?7p~ 



r 

Jo 



dx 



cosh(v / x 2_ +~o 2 ") + cosh(p) 
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(22) 



To obtain the 0(p°) term in the bracket of Eq.(21) we used the integration formula [16] for 
small a 2 : 



dx 



Vx 2 + & 2 (l + e 



lnf^j + 7 + 0(a 2 ) 



(23) 



Integrating Eqs.(21) and (22), at sufficiently high temperature, we obtain 



Vr(v, p) 
N 



1 

7T 
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2 1/1 7C(3) 2 2 

* + 2{-R 2 -n) P + 4^ P + 



(24) 



up to a field-independent constant. Since the 0(5"°) term in the bracket of Eq.(22) is 
exact, Vt{(t, p) has no term higher than the p 2 term that consists of p fields only. At high 
temperature, the Z 2 chiral symmetry is always restored. However, the behavior of the charge 
conjugation symmetry is the same as that at zero temperature. Hence, charge conjugation 
symmetry is not restored at high temperature, for R 2 > R 2 . 

Until now in order to find the effective potential Vr(a, p) we have used the mean field 
approximation (MFA) by taking a and p as constant fields. This corresponds to the leading 
approximation in the 1/N expansion where the a and p loops (i.e., the fluctuations of the 
cr and p fields) are not included. Let us check the validity of our calculations. In the large 
N limit MFA is good, while for large but finite N it may fail due to the contribution from 
kinks in two dimensions [pjjnjfl . 

At first consider the case of the usual two-dimensional GN model (without the random 
chemical potential) where MFA predicts a wrong critical temperature To 0) [|l7]]. For < 
T < T the MFA effective potential is the double well with two degenerate minima at M(T) 
and — M(T), the solutions of the gap equation. Due to quantum tunneling between two 
degenerate minima, the system has kink solutions alternating between M(T) and —M(T). 
They have higher energies than the constant solutions M(T) and —M(T). The Helmholtz 
free energy F is related to the internal energy U and the entropy S by F = U — TS. Since 
F — U at zero temperature, the constant solution is favored and MFA is expected to be 
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valid. If we consider the contribution from kinks explicitly, we find the average (or blocking) 
potential [19j which has a plateau between M(= M(0)) and — M ||17|| . Even in the presence 



of an arbitrarily small external field this potential is tilted to favor M or —M. Hence the 
contribution from kinks does not change the physics materially and MFA is qualitatively 
valid at zero temperature. For < T < T the number of kink configurations is sufficiently 
large to gain enough entropy and so their probability is overwhelming. Since the region of 
a = M(T) will, on the average, have the same weight as those of a = —M(T), we have 
(a) = 0, which indicates the breakdown of MFA. For T > Tq the system is in the symmetric 
phase (M(T) = 0) in MFA and thus has no kink solutions. Consequently, for the two- 
dimensional GN model MFA is good only at zero and high temperature (T > T ) and, by 
the formation of kinks, the true critical temperature turns out to be zero. 

Now let us examine the validity of MFA for the two-dimensional GN model with the 
random chemical potential Eqs.(3) or (4). For the purpose of the present paper we consider 
only the cases of zero and sufficiently high temperature. At zero temperature the MFA 
effective potential V (a, p) has degenerate minima at (|cr|, |p|): (M, 0) for < R 2 < R 2 and 
(0, oo) for R 2 > R 2 (see Fig. 2 (a) and (b)). Hence for < R 2 < R 2 our system may have 
kink solutions for a alternating between M and — M, but no kinks for p. In this case the 
situation is similar to that of the usual GN model in the previous paragraph. Thus it is 
expected that MFA is good at zero temperature and for < R 2 < R 2 . For R 2 > R 2 the 
MFA effective potential is unbounded from below and there is no tunneling between two 
degenerate ground states because of infinitely high barrier. So a and p do not have kink 
solutions. At sufficiently high temperature the a and p fields are decoupled from each other 
in Vr{cr,p) and can be treated separately. By the restoration of the Z2 chiral symmetry 
in MFA, a has no kink solutions, irrespective of the value of the disorder strength. For 
< R 2 < R 2 the charge conjugation symmetry is preserved in MFA and so no kinks for 
p. For R 2 > R 2 the situation is the same as that at zero temperature. For R 2 = R 2 the 
MFA effective potential for p vanishes. As R 2 is tuned through R 2 , the charge conjugation 
symmetry undergoes a first-order phase transition from the symmetric phase ((p) = 0) to 
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the broken phase ((p) = oo or — oo), following positive or negative values of (p) according 
to the value of (p) in the broken phase. So it is reasonable to assume that p has no kink 
solution at R 2 = R 2 . As a result, MFA is reliable at zero and high temperature for all values 
of the disorder strength. 



IV. ORIGIN OF CHARGE CONJUGATION SYMMETRY NONRESTORATION 



In this section we discuss the mechanism of SNR for the charge conjugation symmetry 
conceptually. For convenience set o = in Eqs.(4) and (5) because SNR is the effect of 
the random chemical potential. That is, we neglect the four-fermion interaction (\l/\l/) 2 and 
consider the system of free massless Dirac fermions in the presence of the random chemical 
potential. To leading order in the 1/N expansion the finite-temperature effective potential 
Vt(0, p)/N consists of two parts: (i) a term from probability distribution (p 2 / (2R 2 ) = P#(p)) 
and (ii) the grand free energy (or grand potential) for free massless Dirac fermions at constant 
chemical potential p (= Qt (p))- 

So it is essential to conceptually determine the value of the chemical potential that f^r(p) 
favors. By the symmetry of f^r(p) we can restrict ourselves to the positive chemical potential 
without loss of generality. At first, consider the case of zero temperature. According to 
Fermi-Dirac statistics, fermions fill all the energy levels to the Fermi energy (=chemical 
potential) and antifermions are suppressed. Hence, the larger chemical potential, the larger 
charge (number) density (= fermion number density— antifermion number density). This 
result is retained at nonzero temperature. Since the grand free energy density is minus the 
pressure p0| , it is a decreasing function of the (positive) charge density. Consequently, the 
large charge density (i.e., large chemical potential) is preferred and, for all temperatures, 
^r(p) is minimized at large chemical potential (\p\ — > oo). This implies SNR for the charge 
conjugation symmetry in the open system of free massless Dirac fermions. Moreover, we can 
guess the functional form of f^r(p) by dimensional analysis: At zero temperature, Qo(p) oc 
— \p\ d in d dimensions. At sufficiently high temperature, fir(p) oc —Tp 2 + 0(p 4 /T) in three 
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dimensions and £It(p) oc — p 2 + 0(p 4 /T 2 ) in two dimensions, up to a field-independent 
constant. These qualitative results can be checked explicitly from Eqs.(15), (19), (20) and 
(24). 

Now let us consider the contribution Pr(p) from the probability distribution of the ran- 
dom chemical potential. In three dimensions, for large \p\, flx(p) exceeds Pr(p) at zero and 
sufficiently high temperature, irrespective of the value of the disorder strength. Therefore, 
the initially broken charge conjugation symmetry is not restored at high temperature. In 
two dimensions Pr(p) is comparable to ^t(p)- The probability distribution of the random 
chemical potential for weak disorder (small R 2 ) is dominated at p = and the charge con- 
jugation symmetry is preserved at zero and high temperature. However, for strong disorder 
(large R 2 ) all values of the chemical potential have small probability density and so the fate 
of the charge conjugation symmetry is determined by ^t(p)- Thus, in this case, the initially 
broken charge conjugation symmetry is not restored at high temperature. 

V. CONCLUSIONS 

In the present paper we examined the symmetry behavior of the Gross-Neveu model with 
random chemical potential which is equivalent to the four-fermion model Eq.(3). We used 
the leading approximation in the 1/N expansion (i.e., the mean field approximation). Our 
model has the charge conjugation symmetry as well as the Zi chiral symmetry. The initially 
broken Z 2 chiral symmetry is always restored at high temperature. In three dimensions, 
the charge conjugation symmetry that is broken at zero temperature, is not restored at 
high temperature, irrespective of the value of the disorder strength R 2 . In two dimensions, 
at zero temperature, the charge conjugation symmetry is not broken for weak disorder 
(0 < R 2 < R 2 {= tt/2)), but broken for strong disorder (R 2 > R 2 ). Therefore, our system 
exhibits a quantum phase transition at R 2 = R 2 as the value of R 2 is varied. For any given 
value of R 2 the high-temperature behavior of the charge conjugation symmetry is the same 
as its zero-temperature behavior. Hence the charge conjugation symmetry remains broken at 



12 



high temperature (i.e., symmetry nonrestoration) for R 2 > R 2 . By examining the existence 
of the kink solutions we checked that the mean field approximation is reliable even in two 
dimensions at zero and high temperature. 

In addition, we discussed our results on charge conjugation symmetry nonrestoration 
conceptually, after neglecting the four-fermion interaction (*1>\I/) 2 for convenience because 
symmetry nonrestoration is the effect of the random chemical potential. The behavior of 
the charge conjugation symmetry is determined by the competition of two terms in the 
finite-temperature effective potential: (i) the term p 2 / (2R 2 ) from the Gaussian distribution 
for the chemical potential p that favors the symmetric phase ((p) = 0) and (ii) the grand 
free energy for free massless Dirac fermions which favors the broken phase (\(p)\ — > oo). 

As further work, it would be worthwhile to perform next-to-leading order calculations 
and consider a non-Gaussian distribution for the random chemical potential. 
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FIGURES 

FIG. 1. The zero-temperature effective potential Vo/(NM 3 ) in three dimensions as a function 
of a/M and p/M for R 2 M = 10, in the case of the Z2 chiral symmetry breaking (G 2 M = —it). 

FIG. 2. The zero-temperature effective potential Vq/{NM 2 ) in two dimensions as a function of 
a/M and p/M, for (a) R 2 = 1 and (b) R 2 = 10. 
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FIG. 1 . S. I. Hong and J. B. Kogut 




FIG. 2 (a). S. I. Hong and J. B. Kogut 




FIG. 2 (b). S. I. Hong and J. B. Kogut 



